Critical current in Nb-Cu-Nb junctions with non-ideal interfaces 
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We report on experimental studies of superconductor (Nb) -normal metal (Cu) -superconductor 
(Nb) junctions with dirty interfaces between the different materials. By using a set of simultaneously 
prepared samples, we investigated the thickness dependence as well as the temperature dependence 
of the critical currents in the junctions. Good agreement between the decay of the measured critical 
currents and theoretical calculations was obtained without any fitting parameters. 



a large ratio of interface to normal-metal resistances. We 
demonstrate good quantitative agreement with the the- 
ory of Kupriyanov et al. A in the appropriate limit. 



II. THEORETICAL BACKGROUND 

For the long SNS junction, the one dimensional Us- 
adel equations 8 are used with the appropriate boundary 
conditions. In the zero approximation in the Matsubara 
frequencies, the critical current in a junction with ideal 
interfaces can be expressed as follows 2 : 



Ir. = 



64nk B T 



A 2 (T) 



eRj [7TkBT + n+(2n(TTk B T + Q)) 1 / 2 } 2 ^ N 



e 5 » 
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I. INTRODUCTION 

The non-dissipative current, as well as other super- 
conducting properties, are induced in the normal metal 
by an electron-hole pair that is formed during the An- 
dreev reflection at the superconductor (S) -normal metal 
(N) interface 1 . The finite length over which the electron- 
hole pair remains correlated determines the character- 
istic decay length, £at, of the superconducting current 
in SNS junctions. The normal-metal thickness, L, and 
temperature dependence of the super-current in the long 
SNS junctions, L > £jv, has been thoroughly investi- 
gated over the last few decades, both theoretically 2-4 and 
experimentally 5-7 . 

However, only in the case of junctions with clean in- 
terfaces, was reasonably good agreement between the 
theory 2 and the experiment 7 found. In practice, there is 
no way to produce vanishing boundary resistance. Nev- 
ertheless, the ratio of the latter to the resistance of the 
normal metal in the junction could be made much smaller 
than unity, as was the case for the lateral junctions in Ref. 
7. In this limit, the theories developed for the ideal inter- 
face junction are valid. It should be emphasized, that the 
high quality interface obtained during in situ sequential 
deposition of two materials could result in large values 
of the ratio provided that the resistance of the normal- 
metal in the junction attains a low value. For non-ideal 
junctions, namely, for devices with vertical junctions, the 
ratio of interface to normal-metal resistances is usually 
larger than unity. Therefore, for such devices the the- 
ory of the ideal interface limit 2 is not applicable. The 
influence of boundary scattering on the critical current 
in SNS junctions has been considered theoretically in the where <tat is the normal metal conductivity, and Tb — 
past 3 ' 4 , However, so far, there is no experimental data IklS. Here, 7s is a parameter related to the trans- 



where A(T) is the superconductor energy gap, Q = 
^(irk B T) 2 + A 2 (T)), and £ N = J ' the thermal 
length. 

In the case of non-ideal interfaces (transmission less 
than unity) between the normal metal and superconduc- 
tors, the Eilenberger equations 9 must first be solved for 
the vicinity of the boundary, after which it is possible to 
write the boundary conditions for the Usadel equations. 

Following Ref. 4, the resistivity of the structure in the 
normal state can be expressed as 



(i + 2r B ) 



(2) 



which quantitatively verifies the theoretical prediction. 
The most remarkable difference between the theoretical 
predictions for the ideal and non-ideal junctions is re- 
flected in the dependence of the critical superconducting 
current, Ic, on the total junction resistance in the nor- 
mal state, Rj. For the ideal junctions, where normal- 
metal resistance is the main contributor to the total Rj , 
Ic <x RJ 1 ■ In the opposite limit, where the main con- 
tribution to Rj arises from the interface resistance, the 
supper-current is proportional to RJ 2 - 

This paper is devoted to experimental studies of the 
critical current in vertical Nb-Cu-Nb long junctions with 



parency of the interface. 

In the limit of SN boundary with small transparency 
7b >> (^■) 1 / 2 , the critical current in the first order of 
the Matsubara frequency is given by 4 : 
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where T c is the superconductor critical temperature. 
For the low transparency junctions, namely Tg >> 1, 
the critical current can be expressed in terms of measur- 
able quantities as follows: 



1 



Ir = 
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where pn is the normal metal resistivity and = 
Eq. 4 is valid only for L > (jy. For L > £n 



hD 



27rk B T c ■ 

the contribution of higher harmonics of the Matsubara 
frequencies should be included (see Eq. (34) of Ref. 4). 



III. SAMPLE PREPARATION AND 
MEASUREMENT SET-UP 

The samples were fabricated as long triangular prisms. 
Two SNS junctions were created along the prism in 
a SQUID like geometry as shown schematically in 
Fig. 1 . This design allows the observation of Little-Parks 
oscillations 10 which give two advantages. First, it is pos- 
sible to ensure zero magnetic flux through the junctions. 
Second, accurate measurements of the normal resistance 
and the cross section area of the junctions is possible. 
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FIG. 1. Schematic drawing of a sample. A cross section of 
the structure is given on the right side. 

The SNS junctions were fabricated on top of V-shaped 
substrates, produced by wet-chemical crystallographic 
etching 11 of GaAs planar substrates. 

The V-groove substrate was covered with Nb and Cu 
films, which produced the two bottom side-walls of the 
prism. Nb films were sputtered using a magnetron gun 
and covered in situ with a Cu layer by thermal evap- 
oration, to prevent Nb oxidation. An additional Cu 
layer was e-gun evaporated in a separate vacuum cham- 
ber. These layers were completely covered by photoresist. 
Oxygen plasma etching was applied in order to expose the 
top parts of the V-groove. The subsequent depositions of 
the Cu and Nb films, as described above, completed the 
fabrication of the prism samples. The crossing geometry 



of the top and the bottom films allowed four terminal 
measurement of the junction below the critical tempera- 
ture of the Nb films. 

Nine sets of samples were prepared simultaneously in 
order to assure identical interfaces at all junctions. The 
samples of each set had an identically thick Cu layer. 
Variation of the Cu layer thickness between the different 
sets was obtained by a specially designed shutter, which 
exposed the samples in sequence, so that every set was 
exposed to the evaporating Cu for additional fragments 
of time. Thus, the only difference between the sets was 
the thickness of the Cu layer. We prepared separate Cu 
films simultaneously with the deposition of the Cu in the 
junctions. The measured resistivity of these Cu films at 
4.2K was p N = 5.5 • lO^flm. 

The thickness of each Nb layer was 2000A, whereas the 
Cu thickness varied from L = 2000A to L = lOOOOA. 
The area of the junction, A, , relevant for the critical 
current is determined by the width of the Nb strips, 
WNb — lOyum, and the width of the exposed top part 
of the V-groove, Wt op - The cross-section of the junction 
relevant for the flux dependence of the critical current 
in the junction is given by LWt op - The overall flux de- 
pendence of the critical current in the SQUID sample is 
determined both by the area of the junction and the to- 
tal cross section of the SQUID (Asquid — 6.5pm 2 ), as 
given in the following expression 12 : 



cos ■ 



7T<E> 
$0 



sin(7T$j/$o) 



7T$j/$0 



(5) 



where I Cj is the critical current of one junction, while <& 
and $j are the flux in the SQUID and junction respec- 
tively. 
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FIG. 2. Typical magnetoresistance plot as measured in a 
sample with a layer of 800uACu. 
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The measurements were performed in a 4 He cryostat in 
the range 1.3K-4.2K. The critical current was measured 
by passing a DC current with a small AC modulation 
through the sample. The AC voltage, which appeared 
above the critical DC current, was picked up by a lock-in 
amplifier operated in transformer mode. 
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FIG. 3. Typical variation of the differential resistance of 
a SNS junction. The onset of the peak was defined as the 
critical current of the junction. 



IV. EXPERIMENTAL DATA AND ANALYSIS 

We start the description of our data with the magne- 
toresistance oscillations of the samples. Fig. 2 shows a 
typical magnctoresistance plot for a sample with a Cu 
thickness of L = 8000A The short oscillation period, 
AH = 3.1 Gauss, corresponds to the area 6.5/zm 2 en- 
closed by the SQUID, which is consistent with our de- 
signed geometry of the samples. The beating modula- 
tion with larger period and attenuation of the oscillations 
arises from additional Fraunhofcr oscillations of the criti- 
cal current in each junction (see Eq. 5), with typical area 
LW top - Therefore, the analysis of the magnetoresistance 
oscillations allows us to extract Wtop = 1.2/ira, resulting 
in Aj = 12/im 2 . Thus, magnetoresistance oscillations de- 
termine the geometrical dimensions of our samples, given 
in the previous section in the parenthesis. 

The critical currents in our samples were measured at 
temperatures in the range of 1.3K-4.2K using the setup 
described above. A typical differential resistance of a 
sample vs. DC current is plotted in Fig. 3, which clearly 
indicates the sharp onset of dissipation at / = I c . 

Some of the parameters which appear in Eq. 1 and 
Eq. 4, namely Rj = 19mil,T c = 8.0K and p Cu = 
5.5 • 10 fim, are directly measured in our experiment. 
The rest of the parameters, namely £n, Cat an d A can 
be calculated. The diffusion coefficient, which appears in 



£jv and can be calculated from the measured resis- 
tivity of the Cu, using the relation D = (e 2 N(ep)p)~ 1 , 
where N(ep) = 1.56 • 10 47 J -1 is the density of states at 
the Fermi energy for Cu 13 . From the latter, we obtain 
D = 4.5 ■ 10~ 2 m 2 /sec. The energy gap is calculated us- 
ing the known relation for Nb 12 A = 1.9ksT c — 1.3meV. 
Thus, there are no free parameters entering the theoret- 
ical expressions. 

We should keep in mind that Eq. 3 is only valid for 
large values of 75. The value of 73 calculated from Eq.2 
is about 500, and consequently, Tb varies between ~ 50 
for L = 90001 to - 200 for L = 2000l. These values 
justify the applicability of both Eq.3 and Eq.4. 
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FIG. 4. Universal curve of the critical currents (dimension- 
less units), and the theoretical fits with the parameters de- 
scribed above. The dashed line represents Eq. 4, the solid 
line represents the infinite series (Eq. 34 in Ref. 4) and the 
gray area represents the theory for clean interface junctions 2 
with thickness between 2000 A and 9000 A. 



Since our aim is to compare our experimental data 
with the theory of Kupriyanov et. al 4 , we choose a di- 
mensionless frame of axis (L/£at, I c (T, L)/Iq(T)), where 



Io(T) = 



Wirk B T c 



e N 



These axes are 



eA B?. (7rfe B T) 2 +A 2 (T) £„ 

natural for the I c variation since the leading term in the 
Matsubara frequency expansion (Eq. 4) is just a sim- 
ple exponential dependence, and therefore all our data 
points, presented in this frame of axis, should collapse 
onto a single exponential curve. Fig. 4 shows all our 
data points which quite closely follow (within a factor of 
two) the theory of Kupriyanof et. al. The dashed line in 
Fig. 4 corresponds to Eq. 4, while the solid line accounts 
for all Matsubara frequencies (Eq. 43 in Ref. 4) . 

The consistency with the latter theory is even more 
striking when the data is compared with the theory de- 
veloped for a junction with clean interfaces 2 . Since for 
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this theory, the chosen frame of axis is not natural, our 
data should fall into the gray area. The upper border 
of the area corresponds to L — 2000A while the lower 
border corresponds to L = 9000A. Our data deviates 
by about three orders of magnitude from the theory of 
Zaikin et al. 

It should be noted that there are no free parameters 
in the evaluation of the theoretical formula. The origin 
of the deviation of our data from the theoretical curve 
is not clear, and could arise from several sources. The 
estimate of the area from the Frounhofer oscillations has 
at least ±20% error. The smaller area would increase Iq 
and therefore, improve the fit to the theory. The addi- 
tional source of the error could arise from the the value of 
p in the Cu layer. As explained above, it was determined 
from the resistivity measurement of the Cu film in the 
planar geometry. Although the small grain Cu films usu- 
ally do not show significant anisotropy, one should bear 
in mind that the deposition of our Cu films was made 
with two interaptions (when moving the samples from 
the sputtering chamber to the e-gun deposition cham- 
ber and back). These additional interfaces with partially 
oxidized Cu could make the resistivity of the Cu very 
anisotropic, being much larger in the direction perpen- 
dicular to the interface than the measured values in pla- 
nar geometry. We have no experimental way of verifying 
this. Apart from increasing Iq, and therefore improving 
the fit, the additional interfaces could make some non- 
trivial changes to the boundary conditions of the Usadel 
equations, which were not accounted for by the theory. 

V. CONCLUSIONS 

We demonstrated that the theory 4 of the critical cur- 
rent in a SNS junction with non-ideal interfaces closely 
follows our experimental data. Keeping in mind that we 
have not used any fitting parameters, we find the agree- 
ment quite amazing. We claim that for vertical junctions 
with a clean normal-metal (even for those prepared in 
situ) the theory of Kupriyanov et. al is the relevant one. 
The theory developed for ideal-interface junctions 2 may 



be applicable either for lateral junctions or for vertical 
junctions containing a strongly disordered metal, since 
only in these cases, the interface resistance could become 
negligible relative to the resistance of the normal metal. 
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